arXiv: 1505.0623lvl [cond-mat.mtrl-sci] 22 May 2015 


Direct Method for Calculating Temperature-Dependent Transport Properties 


Yi Liu, 1 Zhe Yuan, 1,2 R. J. H. Wesselink, 1 Anton A. Starikov, 1 Mark van Schilfgaarde, 3 and Paul J. Kelly 1 

1 Faculty of Science and Technology and MESA + Institute for Nanotechnology, 

University of Twente, P. O. Box 217, 7500 AE Enschede, The Netherlands 
2 Institut fur Physik, Johannes Gutenberg-Universitat Mainz, Staudingerweg 7, 55128 Mainz, Germany 
3 Department of Physics, Kinqs Colleqe London, London WC2R 2LS, United Kinqdom 

(Dated: May 26, 2015) 

We show how temperature-induced disorder can be combined in a direct way with first-principles 
scattering theory to study diffusive transport in real materials. Excellent (good) agreement with 
experiment is found for the resistivity of Cu, Pd, Pt (and Fe) when lattice (and spin) disorder are 
calculated from first principles. For Fe, the agreement with experiment is limited by how well the 
magnetization (of itinerant ferromagnets) can be calculated as a function of temperature. By in¬ 
troducing a simple Debye-like model of spin disorder parameterized to reproduce the experimental 
magnetization, the temperature dependence of the average resistivity, the anisotropic magnetoresis¬ 
tance and the spin polarization of a NisoFe 2 o alloy are calculated and found to be in good agreement 
with existing data. Extension of the method to complex, inhomogeneous materials as well as to 
the calculation of other finite-temperature physical properties within the adiabatic approximation 
is straightforward. 

PACS numbers: 72.10.Di, 85.75.-d, 75.47.-m, 72.25.-b 


Introduction .— Measuring the temperature depen¬ 
dence of electrical transport is one of the most important 
and common experimental probes of condensed matter. 
Although a great deal of what determines the tempera¬ 
ture dependence is understood qualitatively [1] , there has 
been virtually no progress in translating this understand¬ 
ing into quantitative, material-specific studies in the past 
twenty years because of the complexity of the theoretical 
formalisms [2, 3]; the lowest order variational approxima¬ 
tion (LOVA) that is the basis for the successful descrip¬ 
tion of the temperature-dependent electrical and thermal 
resistivities of a number of elemental metals [3] has to the 
best of our knowledge not been applied to more complex 
materials. In particular, it has not been extended to the 
study of magnetic materials. The need to be able to do 
so is pressing because current studies of magnetization 
switching involve large threshold current densities that 
are accompanied by substantial Joule heating [4]. 

Inspired by the success of the “direct” ab initio molec¬ 
ular dynamics approach to studying structural and elec¬ 
tronic properties of matter at finite temperatures intro¬ 
duced by Car and Parrinello [5] , we have developed a di¬ 
rect approach to calculate finite-temperature transport 
properties within the adiabatic approximation. For non¬ 
magnetic (NM) materials, we generate “snapshots” of 
a thermally disordered solid [6] and use first-principles 
scattering theory to determine the scattering matrix [7, 8] 
and related properties [9], Fig. 1(a). The results of this 
two-stage procedure are illustrated by comparing the 
calculated and experimentally measured temperature- 
dependent resistivities of the NM metals Cu, Pd and Pt 
in Fig. 1(b). The purpose of this Letter is to underpin 
and extend these extremely promising results by includ¬ 
ing spin-orbit coupling (SOC) [10-14] to determine the 
temperature dependence of the spin-flip diffusion lengths 


l s i for Pd and Pt, and for ferromagnetic (FM) materials, 
to include spin-disorder [11]. 

NM metals .— To describe thermally induced lattice 
disorder in NM metals, we first use density functional 
theory to calculate the dynamical matrix of the bulk 
metal [17-20]. The eigenvalues and eigenvectors that 
result from diagonalizing this matrix correspond to 
phonon energies and vibrational polarizations, respec¬ 
tively. These are used to construct snapshots of corre- 
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FIG. 1. (a) Illustration of the scattering geometry used to 
calculate transport properties. By populating first-principles 
phonon modes, we generate correlated lattice disorder in 
a scattering region (<S) that is connected to semi-infinite, 
crystalline left (£) and right ( IZ ) leads, (b) Temperature- 
dependent electrical resistivities calculated for Cu, Pd and 
Pt. The green dashed lines for Pd and Pt are results ob¬ 
tained without SOC. Experimental data (black stars) [15, 16] 
and the results of LOVA calculations (blue dash-dotted lines) 
[3] are shown for comparison. 











lated thermal lattice disorder by superposing all allowed 
modes in the scattering region populated at a chosen 
temperature [20]. The disordered region is connected to 
ideal leads to perform scattering calculations using the 
Landauer-Biittiker formalism implemented with tight- 
binding muffin-tin orbitals [8] with atomic sphere poten¬ 
tials displaced rigidly with the atoms; see Fig. 1(a). For 
the chosen temperature, the conductance is calculated for 
a number of configurations of disorder and the resistivity 
extracted by varying the length of the scattering region 
as in Refs. [10, 11], whereby the influence of the leads is 
eliminated. This static lattice disorder scheme for elec¬ 
tronic transport is based upon the Born-Oppenheimer 
approximation, which is justified because typical momen¬ 
tum relaxation times for conduction electrons in metals 
(10 -14 -10 -15 s) are so much shorter than the time scale 
of atomic vibrations (10 -12 s). 

Without introducing any adjustable parameters, the 
calculated resistivities of Cu, Pd and Pt, plotted as 
red circles in Fig. 1(b), are seen to be in very good 
agreement with experiment (black stars) [15, 16] indi¬ 
cating that our computational scheme captures the main 
physics of the electron-phonon interaction in NM met¬ 
als. Although SOC has very little effect on the calcu¬ 
lated phonon spectra, including it in the transport cal¬ 
culations for Pd and Pt increases the resistivity and im¬ 
proves the agreement with experiment. In particular, 
above room temperature the calculated resistivity of Pt 
with SOC is 30% higher highlighting its importance for 
5 d transition metals. SOC lifts the degeneracy of energy 
bands and changes the shape of the Fermi surface. By 
allowing spin-flipping, it modifies the phase space of final 
states that can be reached by electron-phonon scattering. 
The SOC strength for d electrons at the Fermi energy in 
Pd (18 mRy) is only about one third of that in Pt (55 
mRy) so that its effect on the resistivity of Pd is much 
weaker. In addition to the SOC parameter for Cu being 
much smaller, the electronic states at its Fermi level have 
mainly s character (l = 0); including SOC has negligible 
effect on the calculated resistivity. For comparison, we 
also plot in Fig. 1(b) the results from calculations that 
used LOVA [1, 21] to solve the Boltzmann equation [3] 
(blue dash-dotted lines). For the past twenty years, this 
has been the state-of-the-art. 

In NM metals like Pd and Pt, the spin-flip diffusion 
length l s f characterizes the relaxation of a longitudinal 
nonequilibrium spin distribution and is an important 
material parameter that enters the description of many 
spin phenomena like the spin Hall effect, spin pumping, 
etc. [13, 22-25]. Its temperature dependence is partic¬ 
ularly interesting since many experiments are performed 
at room temperature; for current-switching experiments, 
Joule heating can increase the temperature of the ma¬ 
terial significantly so electron-phonon scattering is un¬ 
avoidable. By injecting a fully spin-polarized current into 
a disordered NM metal, we are able to obtain l s { using 



FIG. 2. Calculated spin-flip diffusion length of Pd and Pt as 
a function of temperature. Inset: spin-flip diffusion length 
of Pd and Pt replotted as a function of conductivity (solid 
symbols). The solid lines illustrate the linear dependence. 
Experimental values that are either not sensitive to interface 
spin flipping [26, 27] or take it into account [24, 25, 28] are 
plotted for comparison (empty blue symbols). 


an exponential fit to the calculated spin-resolved conduc¬ 
tance as a function of the length of the scattering region 
[10, 13]. The values of l s f we calculate for Pd and Pt are 
shown in Fig. 2. Both decrease with increasing temper¬ 
ature as 1/T. Pd has a larger l s { than Pt at the same 
temperature because its weaker SOC leads to a smaller 
probability of disorder-induced spin-flip scattering. 

Unlike bulk resistivities that are generally well- 
documented, the values of l s { extracted from experiment 
for Pt and Pd exhibit a spread of more than an order 
of magnitude [27-29]. One reason for this spread is the 
neglect of interface spin-flip scattering when interpret¬ 
ing spin-pumping experiments. This leads to a severe 
underestimation of l s f, especially for Pt [13, 24, 25]. An¬ 
other reason is the variable purity of experimental sam¬ 
ples as evidenced by low-temperature resistivities that 
differ substantially [25, 26, 28]. Indeed, it has already 
been pointed out that there is no especially good rea¬ 
son to expect l s f to depend only on temperature-induced 
disorder [25]. Nguyen et al. present experimental evi¬ 
dence for a linear dependence of l s { on the independently 
measured conductivity cr = 1/p [25] in agreement with a 
relationship due to Elliott [30]. We replot l s f as a func¬ 
tion of g in the inset to Fig. 2 and find a perfectly linear 
relation for both Pd and Pt. For comparison, we also 
plot data extracted from experiments that are either not 
sensitive to interface spin flipping [26, 27] or take it into 
account [24, 25, 28]. The good agreement between theory 
and experiment suggests that the Elliott-Yafet mecha¬ 
nism [30, 31] dominates the spin relaxation in NM metals 
like Pd and Pt. 

FM metals. — In magnetic materials, temperature in¬ 
fluences electrical transport by disrupting not only the 
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FIG. 3. Electrical resistivity for bcc Fe calculated as a func¬ 
tion of the temperature. Green triangles: resistivity arising 
from phonon-induced lattice disorder; blue squares: resistiv¬ 
ity in the presence of phonons and magnons that are both 
populated in the scattering region at a given temperature; 
red circles: resistivity with phonon lattice disorder and un¬ 
correlated spin disorder that reproduces the experimental de¬ 
magnetization curve (inset). The experimental values [16] are 
plotted as black stars for comparison. Because of the small¬ 
ness of the anisotropic magnetoresistance in Fe [32, 33], we 
only consider the case of magnetization parallel to the current 
direction. Inset: temperature-dependent magnetization of Fe 
from experiment (black dots) [34] and obtained by populat¬ 
ing magnons (blue dashed line). The red line interpolates the 
experimental values using a cubic spline method. 


translational periodicity of the lattice but also the mag¬ 
netic ordering. Taking bcc FM Fe as an example, we 
carry out the same procedure as for NM metals to ex¬ 
amine the resistivity resulting only from lattice disorder, 
p p h- The results, plotted in Fig. 3 as green triangles, are 
seen to be much smaller than the measured values and 
scale approximately linearly with temperature above 50 
K (as they do for NM metals), whereas the experimen¬ 
tal data show a higher order dependence on temperature. 
This lack of agreement indicates the need to take other 
scattering mechanisms into account. 

The electrical resistivity arising from spin fluctuations 
has recently been studied using various models [11, 35- 
39], none of which is fully ab initio. By analogy with 
the phonon description of lattice disorder just presented, 
we can introduce spin disorder without any adjustable 
parameters by superposing magnon modes that are pop¬ 
ulated as a function of temperature to generate snap¬ 
shots of correlated spin disorder [20]. In the absence of 
a magnon mass, this procedure can be justified by ap¬ 
pealing to the low frequency of thermal spin fluctuations 
compared to typical electron relaxation times. We calcu¬ 
late the spin waves of bulk Fe using the “frozen magnon 
method” [40] and then calculate the resistivity of Fe with 
both lattice and spin disorder, p p h+mg in the same way 
that we calculated p p h- The results are shown in Fig. 3 


as blue squares. Below room temperature, the agreement 
with experiment is much improved. To the best of our 
knowledge this is the first attempt to include both lattice 
and spin disorder in a quantitative, parameter-free study 
of transport in a magnetic material. 

In spite of the improvement, above room temperature 
P p h+mg is lower than the experimental data. To un¬ 
derstand the deviation, we calculate the magnetization 
corresponding to themally occupying the magnon modes 
[40] and plot it in the inset to Fig. 3 (blue dashed line) 
together with the experimentally measured thermal de¬ 
magnetization curve (black dots) [34]. It is seen that 
below room temperature, the calculated magnetization 
reproduces the measured values quite well but gradually 
deviates from the experimental data as the temperature 
increases. As the magnon mode occupancy increases 
with temperature, it gives rise to relatively large cone 
angles, especially for long wavelength modes. The frozen 
magnon method that is based upon a small cone angle 
assumption (analogous to the harmonic approximation) 
is no longer applicable. In addition, magnon modes start 
to interact with each other and deviate from pure bosonic 
character [40]. 

To examine the effect of overestimating M{T) more 
clearly, we switch to the empirical uncorrelated spin dis¬ 
order scheme (disordered moment: dm) introduced in 
[11] to reproduce the experimental magnetization at ev¬ 
ery temperature in addition to the correlated lattice dis¬ 
order obtained by populating phonons. This results in a 
larger resistivity /9 p h+dm (red circles in Fig. 3) that agrees 
remarkably well with experiment. Below room tempera¬ 
ture, the lack of correlation in the spin disorder leads to 
a slightly higher resistivity than what is seen in experi¬ 
ment. We conclude that the underestimation of the fully 
ab initio scheme, p p h+mg, is attributable to the overesti¬ 
mation of the magnetization at higher temperatures. 

FM alloys .— Within the adiabatic approximation, our 
ability to describe from first principles the temperature 
dependence of transport properties is limited by our abil¬ 
ity to characterize temperature-dependent lattice and 
spin disorder. The limitations posed by the harmonic 
approximation of lattice dynamics could be circumvented 
by using ab initio molecular dynamics to generate suit¬ 
able configurations of disorder. In spite of the many 
efforts made to improve the finite-temperature descrip¬ 
tion of magnetism [41-43], this is still an essentially 
open problem. For FM materials with well character¬ 
ized demagnetization curves, we can use the uncorre¬ 
lated spin disorder scheme discussed above. We illustrate 
the effectiveness of this approach by calculating some 
temperature-dependent transport properties of the im¬ 
portant FM alloy Ni 8 oFe 2 o, Permalloy. For simplicity, 
we describe the (uncorrelated Gaussian) lattice disorder 
using a Debye model [11, 44]. 

Using the Debye temperature of 450 K extracted 
from experiment [45] and the experimentally measured 
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FIG. 4. Calculated temperature-dependent properties of the 
fee Ni 8 oFe 20 substitutional alloy, Permalloy, (a) Comparison 
of calculated and measured [50] average resistivities p(T). In¬ 
set: AMR as a function of the temperature. Experimental 
values [51] at 77 K and 293 K are plotted for comparison, (b) 
Spin polarization P = f° r the magnetization 

parallel to the current direction. P becomes only 1% larger 
when the magnetization is perpendicular to the current direc¬ 
tion. The empty (red) circles indicate the spin polarization 
with lattice disorder only. The (green) crosses are values of 
P extracted from current-induced spin-wave Doppler shift ex¬ 
periments on a 20 nm thick film [52] . The black triangles are 
the results of measurements at room temperature for three 
film thicknesses and “bulk” is the value obtained by extrapo¬ 
lation to infinite thickness [53]. 

magnetization [46] to generate configurations of uncor¬ 
related Gaussian lattice and spin disorder, we calcu¬ 
late the temperature-dependent average resistivity p = 
(2 p± + jO||)/3 for Permalloy shown in Fig. 4(a), where 
p± and p|| are the resistivities calculated with the mag¬ 
netization perpendicular and parallel to the current di¬ 
rection, respectively. The anisotropic magnetoresistance 
(AMR) that leads to the spin rectification effect [47] and 
is responsible for the negative domain-wall resistance in 
Permalloy [12, 48] is also calculated and plotted as a func¬ 
tion of temperature in the inset to Fig. 4(a) as the AMR 
ratio (p|| — p±_)/p. Both p and the AMR ratio are virtu¬ 
ally indistinguishable from the best available experimen¬ 
tal data. The AMR decreases from 11% at 100 K to 4% 
at room temperature. This reduction is consistent with a 
picture of SOC-induced scattering, whereby thermal lat¬ 
tice disorder gives rise to weaker anisotropic scattering 
than magnetic “impurity” (Fe in Ni) scattering [49]. 

The degree of spin polarization P of a FM metal, called 
the spin asymmetry in the diffusion theory of spin trans¬ 
port [54] , plays a very important role in many spintronics 
applications [55] but has been quite controversial [56]. It 
has been experimentally measured using different trans¬ 
port techniques, such as magnetoresistance [57], Andreev 
reflection [58, 59], current-induced domain-wall motion 


[60, 61] and spin-wave Doppler shift [52, 62]. However, 
the values reported for Permalloy do not agree with each 
other [52, 53, 57, 59, 62]. Here we focus on the bulk value 
of P = (jf— ji)/{j\+jx) and its temperature dependence 
by examining the projected current densities jf parallel 
to and ji antiparallel to the quantization axis, respec¬ 
tively [63]. Figure 4(b) shows the temperature depen¬ 
dence of P calculated for bulk Permalloy (solid circles) 
for the magnetization parallel to the current direction; it 
increases by about 1% when the current is perpendicu¬ 
lar to the magnetization. As the temperature increases 
from 100 to 500 K, P decreases monotonically from 0.9 
to 0.53. The main contribution to the reduction arises 
from thermally induced spin fluctuations. If we artifi¬ 
cially switch off spin disorder and only include lattice 
disorder, the calculated spin polarization (empty circles) 
is much larger, especially at high temperature. Experi¬ 
ment finds that P depends not only on the temperature 
[52] but also on the thickness of the sample [53], sur¬ 
face scattering apparently depolarizing the current. The 
value labelled “bulk” in Fig. 4(b) was measured at room 
temperature and obtained by extrapolating the results 
found for thin films [53]. The agreement with our bulk 
calculation is excellent [64]. 

Summary. — We have presented a conceptually simple 
“direct” method for calculating temperature dependent 
transport properties based upon the adiabatic approx¬ 
imation that combines first-principles scattering theory 
with temperature-induced disorder modelled in large lat¬ 
eral supercells. The effectiveness of the procedure is il¬ 
lustrated by the very good agreement that we find be¬ 
tween measured and calculated temperature-dependent 
resistivities of Cu, Pd, Pt, Fe and Permalloy. Because our 
scattering formalism includes SOC, its influence on trans¬ 
port in combination with temperature-induced disorder 
can be studied. We find good agreement with available 
experimental results for the spin-flip diffusion lengths of 
Pd and Pt and for the AMR and spin polarization of 
Permalloy. Our calculation of the intrinsic l s f for bulk 
Pd and Pt should help to resolve the ongoing contro¬ 
versy about its value in different situations and provide 
a way to determine l s { from measured resistivities. 

Temperature-induced disorder modeled in this way can 
be used to calculate other properties at finite tempera¬ 
tures such as densities of states, optical excitations, ther¬ 
moelectric effects etc., as long as the adiabatic approx¬ 
imation is applicable. For example, in spin caloritron- 
ics [65], a temperature gradient can be modeled as an 
inhomogeneous phonon and magnon occupation. Our 
methodology makes it possible to study how transport 
changes from ballistic to diffusive with temperature with¬ 
out building any assumption about the nature of the 
transport into the theoretical approach. This can be es¬ 
pecially important for the complex, inhomogeneous lay¬ 
ered structures that are essential for spintronics devices, 
where mean free paths are longer than the “thickness” of 
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When generating thermal lattice (and spin) disorder 
in the scattering region, the atomic (spin) displacments 
should have the periodicity of the lateral supercells. This 
can be done by calculating all phonon (magnon) modes 
with wavevectors q that correspond to reciprocal lattice 
vectors of the real space vectors that describe the super¬ 
cell. 


I. PHONON CALCULATIONS 

To calculate thermal lattice disorder, we use the “force 
constant approach” [SI, S2] . In this method, a (quite 
unrelated) supercell of the material making up the scat¬ 
tering region is constructed, the central atom is displaced 
by a small amount 6 and the forces induced on all of the 
atoms in the supercell by this displacement are calcu¬ 
lated. As long as it is sufficiently small, the forces should 
be linear in S allowing them to be differentiated numer¬ 
ically to form second derivatives of the energy, the force 
constant matrix. By constructing Bloch sums of the force 
constant matrix for arbitrary wave vector q, we obtain 
the dynamical matrix. Because of the strong screening 
in metals, the force field induced by displacing a central 
atom is short ranged. Using elements of the force con¬ 
stant matrix calculated with a 5 x 5 x 5 superccll yields 
well converged phonon dispersion relations. 

In practice, we use the QUANTUM ESPRESSO density 
functional theory code [S3] based on plane-waves and 
pseudopotentials in combination with the Perdew-Burke- 
Ernzerhof form for the generalized gradient approxima¬ 
tion to the exchange-correlation energy [S4] . The experi¬ 
mental lattice constants were used for all metals. The cal¬ 
culated phonon dispersions for Cu, Pd and Pt are shown 
in Fig. SI. The measured phonon spectra [S5] are in¬ 
cluded for comparison and good agreement between ex¬ 
periment and calculation is found in all three cases. The 
Cu and Pd phonons were calculated without spin-orbit 
coupling (SOC). For Pt, we calculated the phonons with 
(red solid lines) and without (blue dashed lines) SOC. 
Because the interatomic forces are mainly determined 
by the Coulomb interaction between electrons and nu¬ 
clei [S6] and the main features of the electronic energy 
bands are not changed by including SOC, the results (are 
seen in the figure to) lie on top of one another. 


We checked that the density functional perturbation 
theory [S7] yields the same phonon modes as the force 
constant approach. The very good agreement between 
the calculated and measured phonon dispersions seen in 
Fig. SI indicates that the harmonic approximation upon 
which both theoretical methods are based captures the 
most important physics. Should it be necessary to gen¬ 
erate correlated lattice disorder including anharmonic ef¬ 
fects as input to a transport calculation, first-principles 
molecular dynamics calculations could be used. 


II. POPULATING PHONONS AND MAGNONS 

Having obtained the phonon energies w sq and polariza¬ 
tion vectors e sq by diagonalizing the dynamical matrix, 
we are able to populate the phonons in a supercell to 
generate a configuration of lattice disorder for a chosen 
temperature. Here s denotes a particular normal mode. 



FIG. SI. Calculated phonon spectra for Cu, Pd and Pt along 
high-symmetry directions in the fee Brillouin zone. The ex¬ 
perimental data (black dots) are shown for comparison [S5[. 
The calculated phonon dispersions of Pt with (solid lines) and 
without SOC (dashed lines) are nearly the same indicating 
that SOC has very little effect on phonon modes. 
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At a finite temperature T, the vibration of an atom l 
about its equilibrium position R; can be described by a 
linear superposition of all occupied normal modes, 

u ,(T,t) = ^==^£ sq A sq (r)e^ R <-^ t+ ^), (SI) 

sq 

where N q is the number of wavevectors q compatible with 
the lateral supercell used for the scattering region and M; 
is the mass of atom l. </> sq is a random phase of the nor¬ 
mal mode w sq ; varying it allows us to generate different 
configurations of thermal disorder. In our frozen thermal 
disorder picture, we can set the time t to be zero without 
loss of generality. The vibrational amplitude A sq (T) is 
determined by occupying the phonon mode at the tem¬ 
perature T according to quantum statistics. Specifically, 
the quantity w^qAgq/2 should equal the total energy con¬ 
tributed by the sq phonon mode. 

Magnon modes for bulk Fe are calculated using the 
“frozen magnon method” introduced by Halilov et al. 
[S8] and populated as a function of temperature to gener¬ 
ate snapshots of correlated spin disorder by analogy with 
the phonon case. At a temperature far below the Curie 
temperature, the occupation of a magnon mode q results 
in a small polar angle 0 q of local magnetic moments with 
respect to the global quantization axis, i.e., 


M s z/,2 \ = »q( T ) 

2 g^B { N q 


(S2) 


Here M s is the saturation magnetization, g is the Lande g 
factor for the electron, taken to be 2, g,B is the Bohr mag¬ 
neton, () denotes thermal averaging and N q is the total 
number of magnon modes. The temperature-dependent 
occupation of the magnon mode n q (T) follows Bose- 
Einstein statistics. The final polar angle of the magnetic 
moment on every atom results from the linear superpo¬ 
sition of $ q for all contributing magnon modes. 


III. NUMERICAL DETAILS 

The Kohn-Sham potentials in the atomic spheres ap¬ 
proximation (ASA) are calculated self-consistently with¬ 
out SOC using the tight-binding linear muffin-tin or¬ 
bital (TB-LMTO) method [S9]. Experimental lattice 
constants are used throughout. For the slab of collinear 
Ni 80 Fe 2 o binary alloy sandwiched between Cu leads, ASA 
potentials for Ni and Fe are calculated without SOC us¬ 
ing the coherent potential approximation [S10, SI 1] com¬ 
bined with a surface Green’s function method [SI 1] which 
is also implemented with TB-LMTOs. In the surface 
Green’s function calculations, the two-dimensional Bril- 
louin zone corresponding to an fee (111) lxl interface 
unit cell is sampled with a 120x120 grid of k points. 

SOC makes a negligible contribution to the self- 
consistent Kohn-Sham potentials and is taken into ac¬ 
count adequately in the transport calculation using a 


Pauli Hamiltonian approach [SI2]. Such a perturbative 
treatment has been successfully applied in first-principles 
calculations of Rashba splitting [SI3], Dzyaloshinskii- 
Moriya interaction [S14], and our own calculations of 
magnetocrystalline anisotropy [SI2], resistivity and mag¬ 
netization dissipation [S15-S17]. For the same reason, 
the magnon dispersion is calculated without SOC be¬ 
cause it is essentially determined by the exchange interac¬ 
tions; the magnetic anisotropy energy of Fe and Ni 8 oFe 2 o 
is tiny and can be safely neglected. 

The scattering matrix is determined using a “wave- 
function matching” scheme [S18] also implemented with 
TB-LMTOs [S19]. For magnetic materials at a finite tem¬ 
perature, the spin-dependent potentials are rotated in 
spin space [S20] so that the local quantization axis of ev¬ 
ery atomic sphere conforms to the required spin disorder. 
The matrix elements of the Pauli Hamiltonian are eval¬ 
uated using the local quantization axis. We performed 
numerical tests with lateral supercell sizes up to 10x10 
and found that good convergence could be achieved using 
5x5 and 4x4 supercells for transport along fee [111] and 
bcc [001] directions, respectively. The two-dimensional 
Brillouin zones of the 5x5 supercell for fee (4x4 for bcc) 
metals are sampled with 32x32 (28x28) k points, which 
are equivalent to 160x160 (112x112) k points in the cor¬ 
responding lxl Brillouin zone. 



FIG. S2. Area resistance of Pt calculated as a function of the 
length of the diffusive Pt (black dots) using a 5 x 5 lateral 
supercell. The disordered region of length L is connected to 
two semi-infinite perfectly crystalline Pt leads and was con¬ 
structed by populating phonon modes using T = 300 K. The 
red bars show the average values and the standard deviation 
from averaging over more than five random configurations at 
every length. The solid blue line is the linear least squares 
fit. The empty green diamonds are resistances calculated by 
integrating the configuration averaged transmission over the 
energy window defined by —df/de where / is the Fermi-Dirac 
distribution function with T = 300 K. The error bars for the 
diamonds are smaller than the symbol size and hence not 
shown. 
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As a typical example, we plot in Fig. S2 the calculated 
area resistance of Pt as a function of the length of the 
disordered region. The empty green diamonds show the 
resistance obtained by configuration averaging the trans¬ 
mission as a function of the energy and then integrating 
over the energy window defined by the derivative of the 
Fermi-Dirac distribution function with T = 300 K (green 
diamonds) rather than T = 0 K (red bars). The Fermi 
smearing has little effect, partly because the lattice disor¬ 
der already smooths the density of states near the Fermi 
level so that the conductance varies slowly over an en¬ 
ergy range of a few For this reason, the scattering 

matrix was only evaluated at the Fermi level in the re¬ 
mainder of this work. 

Fig. S2 exhibits Ohmic behavior, i.e. the resistance is 
proportional to the length of the disordered region and a 
resistivity value of p = 10.0 ±0.3 /if 1 cm is extracted by a 
linear least squares fit. The extraction does not depend 
on the properties of the leads since the Sharvin resis¬ 
tance and other properties of the leads only contribute 
to the intercept of the linear fit. The computing time 
scales linearly with the length of the scattering region 
and quadratically with the size of the lateral supercell. 
Calculating a single configuration of the longest scatter¬ 
ing region shown in Fig. S2 requires about one hour on 



FIG. S3. Resistivity of Fe calculated with lattice and/or spin 
disorder. p p h (up-pointing, green triangles) is calculated with 
lattice disorder only, obtained by populating phonon modes 
while (a) p mg (left-pointing, violet triangles) is calculated with 
spin disorder only, obtained by populating magnon modes. 
Pph+mg (solid blue squares), obtained with both lattice and 
spin disorder simultaneously, is seen to be greater than the 
sum p p h ± p mg (empty magenta squares), (b) Instead of cal¬ 
culating spin disorder by populating the magnon spectra, pd m 
(right-pointing, orange triangles) is calculated with spin disor¬ 
der described using uncorrelated disordered moments. p p h+dm 
(red solid circles) is calculated with lattice disorder described 
in terms of phonons and spin disorder described in terms of 
uncorrelated disordered moments. p p h+dm is greater than the 
sum Pph ± Pdm (empty black circles). 


a supercomputer node with 32 cores and 256 GB mem¬ 
ory; the calculation parallelized perfectly over the two 
dimensional 32x32 k points summation. 


IV. DEVIATION FROM MATTHIESSEN’S RULE 


For ferromagnetic Fe, we can examine the validity 
of Matthiessen’s inequality by comparing the sum of 
the partial resistivities arising from lattice and spin dis¬ 
order separately to the total resistivity obtained with 
both types of disorder present simultaneously. p p h and 
p mg in Fig. S3(a) are the resistivities calculated with 
only phonons and magnons populated in the scatter¬ 
ing region, respectively. The sum p p h + p mg is smaller 
than the resistivity p p h+mg calculated with both phonons 
and magnons present simultaneously. The same con¬ 
clusion can be drawn for the case using uncorrelated 
spin disorder depicted in Fig. S3(b). Specifically, both 
(Pph ± Pmg)/Pph+mg and (p P h ± Pdm)/Pph+dm are about 
0.9 in agreement with a very recent calculation [S21] . 
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